www.gusucode.com > matlab从零到进阶程序与数据 > matlab从零到进阶程序与数据/第4章 MATLAB程序设计/exampleIntLimit1.m
function sol = exampleIntLimit1(a,e,l) %用嵌套函数表示被积表达式 function f = fun1(beta) f = a.*(1-e.^2)./(1-e.^2.*sin(beta).^2).^(3/2); end %调用fzero求满足条件的beta0值 sol = fzero(@(beta0)fun2(beta0,l,@fun1),3); end %用子函数表示积分 function g = fun2(beta0,l,fhdle) g = quadl(fhdle,0,beta0)-l; end